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Abstract 

In the limit of £ ~ a a /a T — > 00 the gluodynamics without the 
magnetic part of action ( Sm ~ 1/?) is considered as a self-contained 
model. The model is studied analytically in the continuum limit on 
an extremely large lattice ( N T — > 00 ). Scaling conditions for critical 
temperature and string tension are considered. The model shows trivial 
( g 2 ~ a T ) asymptotic freedom in the case of continuous gauge groups 
and nontrivial one ( g 2 ~ l/lnl/a T ) for discrete groups. 
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1 Introduction. MC experiment and some prob- 
lems of renormalisation procedure in LGT 

Since the confirmation of the approach to the continuum limit is of partic- 
ular importance in LGT, the systematic scaling analysis at larger (3 = ^p- 
seams inevitable . In lattice QCD the non-perturbative aspects of the the- 
ory are of primary interest, but the renormalisation-group techniques are 
added , as a rule, in a perturbative way. Feasible enough description of the 
running coupling behavior by field theory methods allows to establish a di- 
rect connection with the standard formalism of perturbation theory and the 
renormalisation procedure in lattice gauge theory (LGT). Now the Callan- 
Symanzik beta function 

Pes (9) = = -bog 3 - hg 5 - b 2 gl - ... (1) 



which leads to 



with 



g 2 = b In — |— + ^- In In — |— + .. 

aA L b aA L 



UN , 34 / N . 

^i = tU (3) 



u 3 16tt 2 ' 1 3 V16tt 2 , 

is computed in a perturbative theory with an accuracy even greater than 
that required in LGT at present. Outside the weak coupling region, how- 
ever, the scaling behavior is described by the beta function, which , generally 
speaking, is unreachable for the perturbation theory. 

In the lattice gauge theory, the conditions, which provide asymptotic 
scaling in a strong coupling regime were established long ago (see, e.g.,(l|]). 
In the case of anisotropic lattice a a j a T ^ 1 such conditions can be written 
as 

At,t - Pa,r exp (a a a a:T M% T ) -> /3 CTjT = const; M CTiT = const (4) 



At,t - At,t exp (a a>T m ajT ) -> a<T = const; m a>T = const (5) 

Although and @ guarantee the scaling behavior of mass gaps m ayT 
and string tensions a a ^ T they cannot be satisfied simultaneously. In both 

cases the Callan-Symanzik beta function @cs — Pcr, T hi ^7 turns to zero at 
the finite value of bare coupling, which scarcely can be regarded as a 'strong' 
one. Moreover, since some physical quantities have different strong coupling 
expansions, it is hardly possible to adjust the couplings in such a way, that 
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all of them will remain simultaneously invariant under changing £ Q . All 
this makes us doubt that the continuum limit may be realized in the strong 
coupling area. 

Even being computed in a weak coupling area Eq. (|]) may differ in dif- 
ferent approaches. In particular, only two leading terms in (|l]) are universal. 
Such disturbances are annoying but hardly significant, because the goal of 
lattice calculation is to discover the value of some physical observable when 
the UV cutoff is taken very large, so that the physics is, indeed, in the first 
term of (||) . Everything else is an artifact of calculation. 

However, although the continuum limit (a —* ) in asymptotically free 
theories corresponds to g — > such theories in the continuum limit do not 
become pertubative ||. For example, while calculating R x R Wilson loop 
expectation value in perturbation theory, the correction to the leading term 
typically has the form 

g 2 In R = g 2 In ► oo; r — > oo, (6) 

so that this correction becomes arbitrary large at large distances. 

At short distances, i.e. for energy scales larger than the heavy-quark 
mass, the physics is perturbative and described by ordinary QCD. For mass 
scales much below the heavy-quark mass, the physics is essentially non- 
perturbative because of confinement. Thus, although a trend in approaching 
to the asymptotic scaling in 2-loop approximation is seen at g 2 ~ 1, one 
cannot exclude that the non-perturbative effects may seriously change Eq. 
(|l|) at extremely small g 2 . 

At the moment, lattice simulations are the most popular way of extract- 
ing truly non-perturbative results from quantum field theories. The lattice 
calculation of a number of physical values, e.g. of the glueball mass, string 
tension, ets., at different g may be used to find the beta function, these 
quantities, however, reflect the long distance properties. The exception is 
the Creutz ratio of Wilson loops in which the self mass contributions are 
cancelled, so it is, indeed, a candidate for the study of beta function. Al- 
though many interesting attempts are made to obtain the QCD running 
coupling from the first principles by lattice computations, however, two old 
problems H are still far from an ultimate solution: the ratios of small loops 
are contaminated by lattice artifacts which results in the systematic error 
v and finite size effects (if the correlation length is comparable with the 
lattice size). 

Detailed analysis || (see also |j] ) of the recent developments in renormal- 
isation technique shows that an important uncertainty arises due to the fact 
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that the results of renormalisation procedure depend on the details of the 
lattice regularization. Since so called bare perturbation theory is regarded 
as unreliable, various recipes, based on mean-field theory or resummations 
of tadpole diagrams, have been suggested to deal with this problem |7|, [|]. 
However, different prescriptions give different results and it is in any case 
unclear how errors can be reliably assessed. 

An interesting method to compute renormalization factors that does not 
rely on bare perturbation theory has been proposed in Q. However this 
method has its own problems, the most important being that the momentum 
p should be significantly smaller than 1/a to suppress the lattice effects, 
but not too small since in this case you cannot apply the renormalized 
perturbation theory with desired confidence. 

A method which may provide general solution of the non-perturbative 
renormalization problem through a recursive procedure [10] is suggested in 
[11, 12, [D], 14]. In such a treatment the scale evolution of the renormalized 



parameters and operators can be studied by changing the lattice size L at 
fixed bare parameters. With an accuracy up to lattice effects the running 
couplings on the two lattices are then related through 

g\2L) = a(g 2 (L)), (7) 

where a is an integrated form of the Callan-Symanzik beta function. The 
recursive procedure has been used to compute the running coupling in 



quenched QCD [11, 12] and the running quark masses [13, H]. An im- 
portant point to note is that the renormalization group invariant quark 
masses M are scheme-independent. Such procedure, however, is not totally 
independent of the initial choice of the scaling functions a(g 2 (L)) in ([?]). 

Let us consider some problems of renormalisation in LGT immediately 
revealed by the MC experiment. Comprehensive studies LGT on lattices 
of up to 48 3 x 16 at finite temperature [[y| and up to 48 3 x 56 at zero 



temperature [16| show that (asymptotic)terms violating the scaling are still 
strong at (3 values accessible for todays numerical simulations. These terms 
manifest themselves most strikingly in the famous dip of the "step function" , 
A/9 = Pes (a) — Pes (2a) referring to the deconfinement temperature [15]. 
This pattern is similar in all SU(N) gauge theories. The strongest sensitivity 
was observed in the N T = 4 case, where the numerical calculations close to 
T c are performed with bare couplings f3 ~ 5.6 -j- 5.8, i.e. in the region of 
the strongest variation of the cut-off with (3 (dip in the A^) |L5|]. Besides, 
with the 'dip' in A^ the presence of even a zero in A^ within this region of 
P may be suspected [17|. Numerical studies showed deviations from (|l|) on 
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entering the scaling region where the correlation length q begins to grow. 
It is especially worth to note, that these deviations are always of such a 
pattern that <f increases faster with f3 than is predicted by (|l|), as if the 
theory approaches the critical point and consequently is not asymptotically 
free |7|. 

Recent analyses [18] showed that scaling violations persist in physical 
quantities such as string tension and hadron masses up to /3 = 6.8 |]J. 
Above the deconfining transition point, the matching procedure inevitably 
suffers large errors. MCRG has been performed by several groups on 16 
lattices in the large f3 region up to 7.2 j|, ^T|, |22|, |23|| . However, these 
results were inconclusive and even controversial. Bowler et al. ]p[| ^] found 
sizeable deviations from 2-loop scaling. Recently Hoek has reanalyzed the 
same data and claimed a very slow approach to scaling 23]. Studies on larger 
lattices with better statistics are required to clarify the scaling behavior in 
the high region of (5. 

Unfortunately, with present computer capacity we are still far enough 
from the scaling region and a-dependent terms are large, so we can hardly 
finally decide if the perturbative beta function ([!]) is appropriate for LGT 
at long distance area. However, if for the calculation of some physical value 
(e.g. critical temperature T c ) one applies the beta function computed within 
LGT for another physical value (e.g. string tension a) essentially better 
results may be expected. Moreover, a itself may be used an instrument 
of phenomenological renormalisation procedure. Indeed, scaling behavior 
for \fo /T c is excellent 0], |2^| . MC computations in 3 dimensions of ra- 



tio t/oz/Tc for the SU(2) and SU(3) cases [p4[ , [25] as a matter of fact, 
show only minor deviations from scaling. This, however, can be taken as an 
evidence in favor of the universality of the renormalisation procedure, rather 
than of the applicability of perturbative beta function. 

Unfortunately, on the basis of todays numerical computations, it is dif- 
ficult to anticipate beta function behavior in the limit of a — > 0, taking 
perturbative calculations as a guideQ. Therefore in the near future, at least 
at this point, analytical methods will be still needed. Though one may well 
believe that deviations from scaling will disappear completely with increas- 
ing lattice sizes, a direct check of this suggestion is rather difficult because 
of limited computer capacity. In our previous article [27] we made an at- 
tempt to compute analytically spatial string tension ^/a in Z (2) and Z (3) 
gluodynamics at finite temperature in spherical model approximation. As 

1 Warnings concerning the enhancement of nonperturbative effects as a — > in the 
presence of power divergences have been issued [E6[ 



5 



it was shown, the spatial string tension demonstrated scaling behavior for 
certain parameter area &o ~ 277 — cos W ) • ^he numer ic value of b^ 
with anisotropy parameter £ ~ 1 agree with the standard value bo = j^^z- 
In this article we try to find analytically the running coupling dependence 
on lattice sizes taking some simple cases and compare such a dependence for 
models with continuous and discrete gauge groups. For this purpose we re- 
consider the well known effective action for finite temperature gluodynamics 
[p8| , that helps us to estimate the critical temperature and running coupling 
behavior in continuum limit on 'large' lattices (N T — ► oo). 

The most serious problem appearing in such approach is the calcula- 
tion of the contribution of the magnetic part (Sm) of the action . It was 
studied in [3C] where it had been established that pertubative corrections 
slowed down the critical temperature dependance on N T . However, there 
are reasons to doubt that the magnetic part contribution can be reliably 
computed perturbatively. Below we give some evidence in favor of the sug- 
gestion that in N T — ► oo limit the situation can hardly be totally cured by 
the computation of the finite set of perturbative series terms. 

Being unable to make such calculations nonperturbatively we consider 
an 'electric toy' model, where all terms of order 1/ £ in QCD action are 
omitted. In particular the magnetic part of the action in this model is 
absent. Although such approach can be regarded as too rough and unre- 
alistic, it may be interesting from mathematical point of view, as an ex- 
ample of a model, where renormalisation procedure can be performed ana- 
lytically entirely within the lattice theory. In fact, very similar approxima- 
tion was partially used for the analysis of pure gauge theories in [42] where 

T = > oo was additionally claimed. As it was mentioned above, along 

with a T —* we consider the limit of N T — > oo , therefore in suggested model 
the finite temperatures area also becomes accessible . Some results obtained 
in such toy model may be useful for general considerations of renormalisation 
procedure in LGT. 



2 Effective action for finite temperature Z (N) glu- 
odynamics 

Here we neglect the contribution of space-like plaquettes i.e. the chromo- 
magnetic part of the action. Of course, one may do it in a strong coupling 

2 Asymmertry (or anisitropy) parameter is defined as £ ~ a a /a T , where a T and a a are 
temporal and spatial lattice spacings. 
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approximation, however, as it was repeatedly pointed out, this regime is 



dangerous to work in and as it has been done in [42] (see also [30] and pip , 
we use the anisotropic lattice 



where £ = \ ^- and g 2 = g T 9a- In the weak coupling region g 2 g 2 



g 2 + 0(g*) @, @soe«^- 

As it is well known, the action for Z (N) gluodynamics in such an ap- 
proximation may be written as 

N T -1 

- S « —Se = /3Re X! X! z ° ( x ' r ) Zn ( x ' r + 1) (x + n, r)* 2 n (x, r)* 

t=0 x,r,n 

(9) 

and having imposed Hamiltonian gauge condition 

jV T - 1 

z (x, r) = l-^+^fl (x) ; fl( X )= l[ z (x, r) ; O (x) G Z (AT) , (10) 

T=0 

we get 

- 5b = /3Re £ (x) z n (x, 0) (x + n)* z n (x, N T - 1)*} - S ch (11) 



The action ( |ll|) is a set of one-dimensional chains 

-S c/l = /?ReW £ z n (x,r)z n (x,r-l)*| (12) 
x,n I t=1 J 

attached by the first z n (x, 0) and last z n (x, iV T — 1) links to the plaquettes 
placed at r = 0. Beside such links, these plaquettes also contain Polyakov 
lines f2 (x) . It is evident that from 

e f3Rez = £ cj . {0) z j. £V = £ exp {ij<p k } = N5i; fk = ^ (13) 

3 [z\ k=0 

follows^ 

1 1 

% iP) = Jr IJ ^ eXp {@ ReZ } = Jf eX P ^ COS + i( Pk3) > ( 14 ) 



fc=0 



3 In the case of, e.g., Z (3) it gives: 3 (x) = e * +2 3 e 2 and S ±1 (x) = e "~| 2 
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and having for each link 

exp j/3 (1) Re [z n (x, r) z n (x, r - 1)*] j 
iV-l 

= ^%(/5)-z n (x,r)^ n (x,r-l)^ (15) 

the sum over all variables except z n (x, 0) and z n (x, iV T — 1) we get the 
expression for (N T — 1) -section chain 

2 exp {ft (N T - 1) Re [(N T - 1) z n (x, r) z n (x, r - 1)*]} 

[z„(x,r)] 

= E O 9 )]^ -1 ^ ( x > ^' ( x » t - V* j > ( 16 ) 

j=0 

that differs from one-section chain simply by the substitution {(3) —* 
[3^ {(3)\ Nt ~ 1 . Summing over the remaining chain variables one easily gets 
the familiar expression 

exp{-5 e// } = 9 ( /3 )^^(^|li||) 7VT r2(x)^(x + nr 

« 9 (^ex P |glM) T J2(x)fi(x + nr| (17) 

and since only small j survive for N T —* oo Eq. (|l7|) leads to 

-S e fftt~p (N T ) ]T « (x) n (x + n)* - N t Fq, (18) 

x,n 

where the effective coupling (N T ) is completely defined by the one-dimensional 
chain and related to the bare coupling (3 by the expression 



and the 'free' term^ is 



F =hiQ? (/3) (20) 



4 Only this term survives if we break the chain at t = to with open boundary conditions 
zo (x, io) ^ (x, to + ^Vt) • It easy to see that the first N T — 1 terms of 'high temperature' 
expansion contribute only in Fq. 
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In the wide area (/? . 1) the functions 9^- (/?) « (/3) (inter alia 9j (/?) « 
for /3 << 1) and therefore 

P (N T ) « ^ (21) 

However, in a very important area of /? >> iV the 9j {(5) values signifi- 
cantly differ from Ij (/?) and thus we obtain 

Mfl „ l + 2e-^— ^)cos(f) 

%(/?) " i + 2e -/ 3 ( 1 - cos ^) 

'2vr 



exp 



2 fl_ cos f^U- /3 ( 1 - cos ^)H (22) 



<9 (N T ) « ( ^ « exp f-2iV T (l - «» f^Y) e^ 1 — L . (23) 



and, accordingly, for iV T >> 1 

-j .ex P (. 
In particular 

, exp(-4iV r e- 2/3 } N = 2 

W) w < } ^ -M /y o ( 24 ) 
exp <^ —3N T e p z > N = 3 



so we come to 3d Z(iV)-spin model with an effective coupling given by ( |23| ) . 
Such models undergo phase transition at the points/3 (N T ) = f3 c . 

Therefore for fixed N T for the critical value of the bare coupling (3 we 

get 

const + In ( x _J (t) ) 

^ * /2 \ (25) 

1 - COS (f 



which is similar to the well-known phenomenological rule 34] 

const 

Pc ~ ( 26 ) 

1 " cos ($) 

5 The critical point for ( [l8| ) is defined as /3 (iV T ) = /3 C (AT) where /3 C (N) are numeric 
constants. In particular f3 c (2) w 0.2288 is the critical coupling of d3 Ising model and 
c (3) « 0.5501 d3 is for Potts model. 
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For N » 1 we cannot cut the series in j at (|14|), but we may 
cos(^)«l-j(^) 2 



% 09) 





/3 /2tt> 
[ 2 {NJ 






2^n=0 


L exp< 


[ 2 ^JVy 


2 n 2 


f 



^ /2tt\ 2 /2tt\ 
2 ^ N J > \N J 



% (/?) 



^3 (1(f) 

Taking into account that for the elliptic function 

oo 

#3 (A; ip) = ^2 ex P {—An 2 + in<p\ 

n=— oo 

f¥ ^ f Up + 27rm) 2 1 
= Va S ex P{ 4A 

m=— oo k ^ 

one can write with good accuracy |3(| 

6» 3 (A; Ri exp {7 (A) cos v? + G (A)} , 
where 7 (A) and G (A) are smooth analytical functions of A 

7(A) « 



2exp{-2A}; A >> 1 

21; A<<1 



G(A) 

and we get for /? (iV T ) 

-ln/3(iV r )«- 
which finally gives 

-ln/3(iV r ) rt 

and 



A » 1 
|lnf; A«l 



% ((3) 



-2\N 7 



2 \Nj 



cxp 



ln/3(iV r )« P«N . 
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If we claim that the critical point (N T ) = (3 C corresponds to the phase 
transition, the dependence of (3 (N T ) on N T must be removed at least in the 
continuum limit by a renormalisation procedure which relates iV T (or a T ) to 
the bare coupling (3. The conditions that provide for such an independence 
are discussed later for some simple cases. 

In the previous expressions we could not avoid somewhat annoying 'gen- 
eralities', but we cherish hope that in what follows those may be regarded 
as 'justified'. Finally we want to stress the fact that the suggested approach 
is not limited to small N T , however, the coupling (3{N T ) in an effective ac- 
tion decrease with N T so we must accordingly increase £ to guarantee that 
P(N T ) »p a . 

3 Effective action for finite temperature U (1), SU (2) 
and SU (3) gluodynamics 

Wilson action for SU(N) gauge group may be written in the static gauge 
U (x, r) = Uq (x) as 

a N T — 1 

~ S E = £ E E ReS P {^o ( x ) U " ( x ' r + X ) U o ( x + n ) U n ( x > r)} (35) 

r=0 x„n 

and after the gauge transformation U n (x, r) — > [Uo (x)] _T U n (x, r) [Uq (x + n)] 
we are sure to get | 

= E (|r ReS P ( X ) ^ ( X ' °) ^ ( X + U ) U n ( X > N r ~ 1)} - (X, n)) 
x,n ^ ' 

(36) 

which, as well as (|l8|) , is given by plaquettes, placed on the last links in 
the temporal direction and represents O (x)-loop interactions with the one- 
dimensional chains 

- s ch (x, n) = ^ReSp j £j U n (x, r) £/£ (x, r - 1) | (37) 

attached to such plaquettes. Having written 

exp {-s ch (x, n)} = d J^j (@) Xj \U n (x, r) f/t ( X) r - 1)J , (38) 

j 

6 The same result may be obtained by imposing the Hamiltonian gauge condition 
U (x, r) = 1 - 5° + S°n (x) with Q (x) = U (x) Nt , from the begining. 
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where 

% 09) = J exp {^ReSp {[/„}} Xj (U n ) dU n (39) 

and Xj (U n ) is the characters of j irreducible representations. So using the 
character orthogonality 

J Xj (U n (x, r) U n (x, r") 1 ) Xk (u n (x, r") C/ n (x, t') ] ) dU n (x, r") 

= ^(l7 B (x,r)tr n (x,70 t ) (40) 

and integrating over all U n (x, r) variables except U n (x, 0) and U n (x, iV T — 1) 
one can obtain for (N T — 1) -section chains 

y exp | ^ p ReSp [u n (x, r) Xj\ (x, r - 1) } | dC7 

= ^d J ^(/3) iV ^ 1 x,(^n(x,0)C/ n (x,iV T -l) t ), (41) 
i 

which, as well as ( |l6|) , differs from a one-section chain simply by (fi) — > 
Substituting([4l|) into ( |36[) and taking into account that 

Xi (fi(x)) Xj (fi(x)t) 



{fi (x) f/ n (x, 0) fi f (x + n) U n (x, 0) f } (if/ 



d i 



(42) 

we can integrate over U n (x, iV T — 1) and finally get the familiar result p8|, p9|] 



exp {-S E } = £ ^P {-iV T (A^ (0) + ln<3? (/?)) } Xj (fi (x)) Xj (fi (x) t 

(43) 

For the gauge groups U(l) and SU(2) we shall have % (/3) = (/?) 
and (/?) = I-y+i (/9) correspondingly. Having made allowance for (|p5[7.13.2(8)) 



7 In (1 + l)-dimensional case we can integrate over f2 (x) and wright for the partition 
function Z = . Qj {f3) NaN ' r . This result is exact because in this case the action does 
not contain the magnetic part. 
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one can write 



T^e^expj-An 2 }; (3 & 1; 



exp 



4n 18 



13.1. 



(45) 



Therefore in the strong coupling region {[3 . 1), we get the result which 
is very similar to that obtained for Z{N) 



4 ] IP) 



InM? vcLln 



2iV 

T 



(46) 



Moreover, taking into account the results |3?| obtained for / \*j (U) e^ x dU, 
Eq. ( |46| ) can be easily generalized for SU(N) case. Therefore, we may con- 
clude that in the region (3 . 1 the effective coupling (3 (N T ) depends on (3 and 
iV T as 



\ const 

For the weak coupling region we can write 



Ag } G9) « A^a 



c 2 (; ) + o(0 



(47) 



(48) 



where C2 (j')is the quadratic Casimir operator, cZ^ = A^ 2 — 1 for SU(N) and 
cIa = 1 for £7(1). The function Ag for £7(1) is given by 



Ai 



hi 



I (/3)J 2/3 V2/37 



+ 



(49) 



and that for SU(2) by 



-3\ _ 9 



(50) 



Similar, but less accurate result can be obtained for SC/ (3) gauge group 
(with j = {Zi;^})- If we apply the asymptotic formula fl45|) to the well 



known expression |37| 

J exp|-/?||xh« 2 ^ = J2 detI ij-j+i-r 



(51) 
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this will give 



J exp|--ReSp(J7)|x W2 d/x 



2 N 



det I h _ J+i ( f) e-^+< 2 ) 2 £ exp -8A 3 f n + ^ ) "I , (52) 



3 

with I2 > h > and 

^ = p( 1 + p)+°^ 3 >- < 53 > 

So, taking into account 

£ exp I - 8 A 3 (n - ^) 2 } = #3 (8A3; 0) + O (exp {-^V) , (54) 

one can easily find 

y exp|-^ReSp([/)|x« 1 /2 d / i 
« ^detVj+^nf^) « exp {-8A3 -CaGi.Za)}, (55) 



,3, 

where C2 (Zj; Z2) = | (Z? + Zf — Z1Z2) — | is the quadratic Casimir operator for 
SU(3) group and we immediately get an expression AgjLgs (/?) that formally 
coincides with (fi"8|) . Therefore, the effective action is given by 



exp {-S E } « [9f E e- 8 ^ ^** (0 (x)) Xj (n (x)t) . (56) 



In the first order in -p the expression (|56|) coincides with that obtained 
in Hamiltonian approach developed in It is worth to note that in 

p8|, |39|] the Hamiltonian gauge Uq (x, r) = 1 can be imposed on all r-links 



in distinction from (36)) because there are no periodic border conditions. 



As a rule such condition is demanded within finite temperature theory |4 



so at this point the approach of |38| , 39] resembles zero temperature theory 



However, as it is known from field theory [El]], if we fix the gauge in a such 
way we would destroy the mechanism that guarantees the Gauss condition 
Qb = (Q b = [5 bc d n - f abc A^j E c n for QFT and Q b = Y? n =-d K for LGT ) 
and a projection operator should be introduced to restore it. 
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In field theory, as well as in LGT [39] the integration variables <pQ (x) in 
the projection operator correspond to the parameters of static gauge trans- 
formations with the generators Q b [f4l|| and finally all looks as if the static 
gauge condition Uo(~k,t) = Uo (x) = exp ji^Q (x) T fe | had been imposed 
from the beginning (as in fl35|)) instead of Uo(x,t) = 1. Therefore, the 
auxiliary variables YI^Iq Uq (x, t) = Uo {x) Nt = exp ^iN T <j)\ (x) T 6 | in the 
projection operator formally correspond to Polyakov loops Q (x) in (|5(|) and 
the lattice length in the temporal direction N T a T corresponds to the in- 
verse temperature. So, at least in this particular case, the border conditions 
appear not to play a significant role. 

Thus, we can use the results |5(j] obtained in Hamiltonian approach to 
sum over irreducible representation j in (|43| ) which, e.g. for SU (2) group, 
gives 



(2j+l)y 
2 



sm 



(2j+iy 

2 



sin % sin ^ 



sm 

^ e -i(j+i)A 

r^-^ 2 ] 2 \ (57) 

f sm T 



and with good accuracy one can write 



2 sin ^ sin 



« exp < 7 cos — cos — + 5 sin — sin — + G > , (58) 

where 7 y^) ,G ^jand <5 ) are smooth analytical functions of A which 

are very similar to given in (|3~I|), and S ( j ^ , G (^) « 7 f £j . Then we can 
finally write 

exp {-S E } « exp {/3 (iV T ) Xj (O (x)) Xi (« (x) f ) + const] , (59) 



with 

In similar manner one may show that in the case of SU (3) gauge group 
effective coupling dependence on (3 and iV T resembles (|60| ) and can be written 
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as 

(N T ) « / (X G N T ) (61) 



with A G given by (gg) , © and (| 

We may conclude that in given approximation for continuous abelian 
U (1) and nonabelian groups SU{2) and SU (3) the effective couplings (N T ) 
are also totally determined by one-dimensional chains, however, contrary to 
Z(N) case (N T ) is changing essentially slower with the bare coupling 
0. 

In Appendices A and B we made an attempt to perturbatively estimate 
the contribution of the magnetic part of the action. First order correction 
in magnetic coupling a turns into zero for symmetry reasons. The most 
important (from our viewpoint) part of the second order correction /?^2 
leads to simple renormalisation of the effective coupling and at least in a 
limit N T — > oo does not change the results drastically. 

More complicated (and possibly more accurate) expression for 0^2 was 
suggested in |30| . As it was established in this paper the additional term 
0^2 makes the movement of T c with increasing N T slower [3C]. Unfortu- 
nately, it is not so easy to anticipate whether it can totally terminate the 
dependence of T c on N T , especially on extremely large lattices, because the 
perturbative expansion in a is hardly applicable outside of the narrow area 



Moreover, if following [30] we assume that second order corrections in a 
radically change the predictions obtained from zero order, the convergence of 
series in a should be reconsidered and one must prove that the corrections 
of orders higher than the second one are insignificant. Anyway perturbative 
series can hardly be regarded as a reliable instrument for the computation 
of effective action in the area a & 1 , N T — > oo (see also footnote 4) 

As it has been already pointed out, a = p| can be made negligibly 
small in comparison with T = p| by changing the anisotropy parameter £. 
Precise analysis of modifications of /3 TjCr with changing £ has been made by 
F. Karsch in Q: 

1 „ { + 0.39832 



T = i(l- 0.27192) + -; a = ; l = V0^. (62) 

2 5 

However, this going into details cannot change the situation drastically for 
large £ (Hamiltonian limit) and large (continuum limit). In this particular 
area one can hardly hope that the second order corrections in a will play 
such a significant role. 
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4 'Electric toy' model 



Since the pioneering work of Svetitsky and Yaffe it is understood that 
all the relevant properties of the phase structure can be encoded in a suitable 
effective action for the order parameter, e.g. the Polyakov loop. The simplest 
and most popular approach to such an action was suggested in |28|, ^] . We 
use the anisotropic lattice (a T « a a ), that helps to avoid the strong 
coupling approximation (one of the basic assumption in 28, ^]).Then we 
try to trace a T -dependance of the effective coupling and compute Callan- 
Symanzik beta function for iV T — > oo. 

As it is pointed out in [^] at high temperatures Tgy f3 a v 1/T$y 
and j3 T v Tgy (hi our denotation @ a v l/£ and j3 T v £). Consequently, 
the functional integral is highly peaked at the configurations in which the 
spatial link variables are static up to gauge transformations and the system 
is described by zero temperature (i-dimensional gauge theory with coupling 
g 2 TsY- The dynamics of the spatial gauge fields is regarded [42| as being 
qualitatively the same at high and low temperatures. 

The naive interpretation of this suggestion might be as follows: the mag- 
netic part of the action Sm, proportional to fi a v l/£ do not play an essential 
role at high temperatures and the effective action, obtained by integration 
over spatial degrees of freedom should not be influenced by the magnetic 
part in the vicinity of the Hamiltonian limit (£ — > oo). However, there are 
serious reasons to believe that the magnetic part of the action may be of 



crucial importance for creating the confining forces |45L 46 , H7l] , even at high 



temperatures [48]. 

Our computation reveals no drastic changes resulting from lower correc- 
tions in p a (see Appendices A and B), however, one would hardly suggest 
with certainty that the picture will be the same for the whole series in j3 a . 
A simple example would prove that such suggestion may be erroneous. In- 
deed, let us consider d2 - Ising model on the anisotropic lattice. Free energy 
density may be written as 

F ( 7ct ; 7t ) = - J^- = In (1 = ^ (1 = ^ + f ^% In / (63) 

where 

/ = (l + -£) (l + 7 2) - 2 7t (l - 7 2) COS if T - 2 la (l - 7 2) COS tp*. (64) 



and 



7(JT = tanh/3 CT)T . (65) 



^In |42| temperature T$y is denned as T$y = \J ir/ia/N T a T 
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It is evident that for any arbitrary small j3 a there is a critical value of (3 T 

(3 C T = -ilntanh/3 CT (66) 

at which the system undergoes a phase transition. Although perturbative 
expansion in (3 a gives a reasonable value of F, the finite series in (3 a do not 
reproduce the singularity of F at (3 T = /3°. In such simple case, the reason of 
the above can be easily traced. The main contribution into integrand at(|63|) 
gives the area of ip a>T v , so 

In / = In (1 + 7r - 2 7r cos ip T ) + v la - )plv 2 + O ( 7 ^) (67) 

with 

1 + 7r + 2 7r cos if T - 2 (1 - 7 r ) cos</v 1 + 5 7r 
■y = v (68) 

1 + 7r — 2 7r COS if T 1 — 7 r 

hence, in the area of 7r ~ 1 the expansion parameter of free energy F is, 
indeed, rj- 7(T rather then 7(J . 

In the considered example, the restoration of real physical picture with 
extending perturbative series in (5 a is very slow. For this reason, one should 
be careful when cutting the series in f3 a in gluodynamics; at some stage 
the magnetic part contribution may start dominating and might completely 
change the final result. Therefore, QCD without the magnetic part should 
be considered as a specific model. 

An effective action (^) of such electric toy model differs from that ob- 



tained in [28, |29| only by initial assumption: we do not use strong coupling 
approximation and work on the anisotropic lattice. An additional assump- 
tion that the magnetic part may be discarded in the limit £ — > 00 , give us 
the right to study weak coupling region and consider the limit a T — > in 
such model. It is evident, that at the same time we can demand a a — > and 
a T N T = T _1 = const < 00. 

Although (|43|) was computed exactly along the line of Svetitsky-Yaffe 
program, discarded magnetic part of the action may lead to serious distor- 
tions. Therefore, it is interesting to compare ( f43| ) with the analytical results 
obtained without such approximation. In |^2[ the analytical solution was 
found for (2 + l)-dimensional U(l) -gluodynamics and a relation was estab- 
lished between (2+l)-dimensional [/(l)-gluodynamics and 2-d Coulomb Gas 
model (CGM), which undergoes a Berezinskii-Kosterlitz-Thoulessf phase 
transition at the point of 2tt 2 ^ = (3° GM « |. 

The model action ( ]43| ) connects (2 + l)-dimensional U(l) -gluodynamics 
with 2d-XY model. As it is known, the 2D X-Y model can be mapped to 
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the 2D Coulomb gas (for recent review see |49|]). Therefore, in exact [}Eq ] 
and approximate approach one comes to similar description of the phase 
structure (2 + 1)- dimensional [/(l)-gluodynamics. The only thing which 
remains is to compare the prediction of [42] for the critical coupling (3 C ~ 
-7^ k, L2 j? 2 (in a T = lunits ) with that one which follows for N T >> j3 

from (|43|), that is with j3 c — (in the same units), where /3^ KT is the 

critical coupling of Berezinskii-Kosterlitz-Thoulessf phase transition. Mote 
Carlo experiment [50| gives fi® KT ~ 1.11. For the toy model, we cannot 
expect more. 



5 Effective coupling and critical temperature 

The effective action (^) for Z (N) - gluodynamics has the critical point at 
f3 (N r ) = f3 c (N T ) where (5 C (N T ) are numeric constants and the critical point 
will depend only on N and the space dimension. If we assume that this 
critical point does correspond to the temperature phase transition, then as 
it follows from (|2l|) this will happen at the fixed temperature T c (in the 
strong coupling regime) when 

(3 \ N - 

= const, (69) 



const 



then by the appropriate choice of constants one can obtain scaling having 
imposed the condition (|5|), but as for the string tension the scaling is hardly 
possible under (||) condition. 

In the weak coupling area, in accordance with (^3|) , and (|4|) we may 

put 



(5{N T ) 



exp{-4iV T e- 2/3 } ; N = 2; 

exp{-3iV T e-^I}; N = 3; (70 ) 



^ exp(-AT T (f) exp{-f(f)' 
so as it follows from 



, 2ir\ 13 (2tt\ 2 



N » 1. 



/2vr\ 2 M^ 2 ) 
N T [ — ) (Aa T ) w ~ — In /3 C ~ const, (71) 



N 
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the critical point (e.g. for Af » 1) will escape undesirable dependance on 
iV T if we demand 

2^"^ 2 " 6oln A^ ; A = const ( 72 ) 

where 

^0 « ^ « - 025iV > ( 73 ) 

which for £ = 1 and iV gives good agreement with bo = j^^i ~ 0.023iV . 
Therefore, we may conclude that 

6 (47T 2 ) 

T c ~ const x a T w ~ const x a~ 1+ ^. (74) 
As it has been shown [27] 

° * exp { _/? (' " cos ^) } " exp r^ 6 ° ln i } (75) 



and ( |74[ ) at given approximation the scaling will take place or, in other 
words, we shall have 

^— = (a T ) ^ \? a) ■ const. (76) 

J c 

In principle, by changing £ one can tune the parameter bo to make T c in 
(|74j) or -^/cr in ( |75| ) tend to the finite value in the continuum limit, however, 
at least, under given approximation it cannot be done simultaneously and 
the ratio ^ will be independent from lattice spacing only for £ = 1. 

6 Continuous groups 

In this case the requirement that the position of critical point (5 (N T ) = j3 c 
corresponding to the temperature T = T c of phase transition should be 
independent of lattice sizes, (at least for iV T —* oo) leads to the condition 

iV r A G -> const. (77) 

For large enough (3 (small a T ) X z g ~ ^j, so (J?7|) means trivial asymptotic 
freedom 

1 g 2 

— — lim — = T/v = const. (78) 
2N a T ^o a T 
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For example, in the case of SU(2) gauge group from ([77]) and (|6(]) we obtain 
the condition 

/3W-7fe 2 )-7f^W c (79) 



N J \N T c 



and consequently 

T N = TcNT 1 (&) , (80) 

where 7 _1 (x) = y is the function inverse to 7 (y) = x, so 7 _1 ^Pcj is simply 
a numeric constant. 

Now we may write for the effective coupling 

(K) » 7 (V (/9c) (81) 

with 

For the calculation of temporal string tension a, higher orders in a T 
are needed in the expansion of g 2 . Indeed, to find a we may compute the 
correlation function (xoXr) between two probes[] 

W=/ exp <!" fi} 7 = ^, (83) 

•/ s - 7 ZnCQSQ'n \2vr/ 2aT 

That can be done in spherical model approximation in the same way as in 
p6[ . The saddle point so is defined from the condition (xoXr) = 1 an d is 
equal to 

so ~ { 3 7c + 8vr 2 7c ( 7c - 7 ) 2 ; 7 . 7c . (84) 
For i? >> 1 one may write 

(XoXr) ~ (27TJ2)- i exp {-a!2} (85) 
and for string tension a we get 

a ~ (7 C - 7) • (7 C - 7) (86) 

or 

a « ^ {T) 2 ~ 2 t (TC) ° (Tc " T) ■ (8?) 
9 Distance R is measured in lattice units. 
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If we claim the independence of a from lattice sizes then we should demand: 

g 2 (T) « { -4Taln (exp _ « a ) ; T<T C . (88) 

This expression guarantees scaling for a in given approximation. How- 
ever, if in continuum limit a — > const then, at least, the second order ex- 
pansion g 2 T in a is needed and, consequently, in this case we must compute 
Pes (<? 2 ) m the higher orders in g. An example of a such computation is 
given in the Appendix C. 



As it is seen from Eq. (117), 0os (y) ~ ~~ § wnen 9 ^ 0> which 



strongly contradict to standard expression (Q) , where 0c s {g 2 ) ~ —bog 3 in 
corresponding area. It looks so as if for extremely small a (and consequently 
extremely small g 2 (a) ) the 'recursive' coupling g 2 (Aa) = F (<? 2 (a); A) might 
be presented as 

g 2 (\a)^c {\)+c(\)-(g 2 {a)) a + ... (89) 

where a is not obligatory an integer, but is independent of A. Asymptotic 
freedom leads to cq (A) = 0, so we immediately get 

0Cs(9) = a^^^rg + O(g 3 ). (90) 

Although the phase structure of the system defined by the action ( |35| ) and 
(|43|) essentially depends upon space dimension d, the effective coupling is 
not sensitive to it and is defined by one-dimensional chain. Therefore, it 
is not too surprising that the toy model shows trivial asymptotic freedom 
inherent to Schwinger model)^]. 

7 Conclusions 

Although in a strong coupling region non-universal behavior of Pes (g) is 
observed, in a weak coupling limit we find the dependence of g 2 on lattice 
spacing for Z(2) and Z(3) gauge groups. This dependence is very close to 
that obtained in standard renormalisation theory. However, the dependence 
on lattice anisotropy parameter £ is too strong compared with [44| and differs 



for spatial string tension and critical temperature. Therefore, either an 
additional procedure must be added to remove the remaining ^-dependence 
in physical values or one must work only at £ = 1. In the last case our results 
should be reconsidered, because we worked far from the area of £ ~ 1. 

10 Another cause of similarity is the absence of magnetic field in Schwinger model. 
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Main approximation, used to obtain an effective action for continuous 
groups is that in a Hamiltonian limit (£ >> 1) the magnetic part of the 
action, (proportional to l/£) can be neglected compared to the electric one 
(proportional to £). Perturbative estimation of the magnetic part does not 
change the results considerably. Therefore, in the limit of N T — ► oo and 
a T —* the toy model shows the trivial asymptotic freedom (g 2 ~ a T ) . This 
strongly contradict the results obtained in the perturbation theory. 

Another nameworthy point is to be underscored: a considerable differ- 
ence between Z(N) and continuous groups claimed by equations (72) and 
( |117j ) . It looks somewhat ridiculous so we would like to discuss it in more 
detail. We argue that such difference between Z(N) and U(l) will not dis- 
appear for any large (but finite) N and will be washed out only at iV — > oo 
, because on their way to continuum limit (j3 — > oo) Z(N) groups inevitably 
pass the point of f3 v ^ 2 = const. With further decreasing a the bare cou- 
pling j3 >> and, consequently, effective coupling (3 (N T ) exponentially 

decrease with (3, as it can be seen from (24) which finally leads to (117) , 
so the dependence of g 2 on a is similar to standard and leads to nontrivial 
asymptotic freedom. In the case of U(l) gauge group (which corresponds 
to N = oo) the point of (3 v ^ ^ is evidently unreachable. Therefore, we 
get different results if the limit N — > oo is taken before the computation of 
partition function, and after that, in other words, the limits N — > oo and 
A^ T — > oo do not 'commutate'. 

Finally we want to stress again that despite the gluodynamics without 
a magnetic part is very far from reality, it may be interesting not only for 
simplicity reasons, but also as an example of a model where the renormal- 
isation procedure may be fulfilled analytically and exclusively within the 
lattice gauge theory. At the same time, it is a good laboratory to study the 
nature of some non-perturbative phenomena in LGT revealing the difference 
between the models with discrete and continuous gauge groups and it is not 
unlikely that this difference will be preserved in more realistic models. 
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_ 2N 

Luructtiuii series 111 \j a 

Z(2) gauge group. 



8 Appendix A. Perturbation series in j3 a = p| for 



Now, for simplicity, we shall confine ourself to the case of Z(2) gauge group 
and, in fact, use the series in j a = tanh/? CT : 

oo 

z(n) = z (n)-J2^n(n) (91) 

ra=0 

The zero order term Zq (O)has been already considered. It is easy to see 
that the first order term Hi is equal to zero. Only the pairs of plaquettes 
having equal spatial coordinates and orientations, but positioned at different 
temporal points r and r + A along the temporal axis contribute into the 
second order E2 . Each of the four chains 1 + O (x) Q (x + n) ^ Nt which cross 
pair links z n (r, x) and z n (t + A,x) of these plaquettes is converted into 

7 A + n (x) n (x + n) 7 ^- A . 

If we denote 

O(x) = Qx; 0(x + n) = 2 ; 
0(x + n + m) = ft 3 ; Q (x + m) = 4 (92) 



and 



^1^2 + 2 3 + O3O4 + ^4^1 
E 4 5 

2 1 + O1O3 + 2 4 + O1O2O3O4 



(93) 



I 

4 

and take into account that I = I 3 and 
1 



IJ (P + qn n n n+1 ) = (f - q 2 ) + Apq (f + q 2 ) l+Sp 2 q 2 l 2 (94) 



. 2 

(p + qn n Q n+1 ) = [p z - q z ^ 

n=l 

after bulky computations we get pT[ : 

jV T — 1 4 



A=l n=l i + 



(95) 

where are the functions of 7 and N T . If e.g. for Z(2) we denote 



7 



exp{-2e- 2/3 } = exp{-2(a T A) 85 °} (96) 
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then for 860 ~ 1 it gives 7^ ~ e £ with e = %rr and for iV T >> 1 we obtain 
sinh 2s — 2e 



Qo 



4e sinh 2 e 




_ (cosh 2e + 5) sinh e — de cosh e 

~ : — — ri : 

4esmh e 

(2ecoshe — 3 sinh e) coshe + e 
8e (sinh e) 4 

The term I contains only the nearest neighbor interactions of n and 
leads to a simple shift of effective coupling. Although the term I 2 also 
includes interactions, that are absent in the standard Ising model, such 



modifications has been extensively studied and it was shown (see e.g. [51] ) 
that they might change the phase structure so drastically that this theory 
would not belong any more to the same universality class as the ordinary one. 
However, if the corresponding coupling (N T ) —ft (N T )is less than (3 (N T ) /6 
(as it certainly is in our case) such interactions does not bring sufficient 
changes into the critical behavior of partition function [51[] . Taking into 
account that only in the 1/8 part of the configurations (ft n = — f2 n +i) > I 
differs from i 2 we put Q + 2Q x l + Q 2 P ~ Qo + (2Qj + Q 2 ) I • It is easy 
to check that in such approximation K + el ~ K exp (j^ ■ so after the 
inclusion of the second order correction we finally get 



with 



ZoiP) v « — y y i + 7^Qo 

« exp [(5(2) (N T ) i + constj (100) 

P {Nt) = 2Ql + 92 (101) 
which evidently leads to a simple shift in effective coupling. 

9 Appendix B. Perturbation series in (3 a for U (1) 
- gauge group 

The zero order term Zq (Q) of the expansion 

Z (Q) = Z (Q) ■ (l + /3 ff 3i (n) + (O) + ... j (102) 
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corresponds to the case (3 a = already considered. It is easy to see that 
the first order term Hi is equal to zero. Computation of the second order 
term: E2 is very similar to that for discrete groups but technical diffi- 
culties increase enormously. The procedure suggested here differs from that 
elaborated in [^] chiefly in technicalities. We are forced to partly sacrifice 
accuracy to obtain a result which appears to us rather simple and transpar- 
ent. 

Only the pairs of plaquettes having equal spatial coordinates and orien- 
tations contribute into the second order term: H2 but they are positioned 
at different points r and r + A along the temporal axis. 

Such additional plaquettes convert one-dimensional chains 

[I j (P)] Nt etffo*-**+»\ (103) 

3 

which cross the pair links U n (r, x) and U n (r + A,x) of these plaquettes 
into 

Y [i j+1 m A ft m Nr ~ A e ^-^\ (104) 



so we get 

l^jlrs { J^yl J cos Wj,r,l,s) 

where 



(105) 



G rr/s = 2cos(^, M ) E p^iVA VW]^ (106) 

A=0 V / \ i / 

and 

'/'iris = j (<^x - V^x+n) + r (l^x+n ~ ^x+n+m) + 
+n+m ^x+m 

) + s (<^ x+m - <^ x ) . (107) 

The term E2 may be easily computed for & & 1 (e.g., A^ T ~ /3 ~ 
In ^- ; a r — ► 0) . In such area, taking into account (jig) , leading terms may 
be written as 

_ 2JV T 

Ooooo = e_i_iu = 2 1 " 6 J (108) 
1-e P 
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The next terms 



©(l) = 2 (0-1000 + @o-ioo + @ooio + 



0001, 



2 — V cos (<p k - <p k+1 ) +Oie 3 (109) 

1 — p 8 u V / 



1-e k 

are of order e . Finalizing we may write 



©0000 + 0(1)- — —ll-e P (l - e p) Vcos((^ fc - p fc+1 ) ) (110) 

1-e 3 V V k / 



and 



I ' u W 

+ T ^2 



^1 + ^^^4 j ex P ^ _ ^(2) cos ~ ^fc+i)^j + O (e 2 & 



(111) 

Therefore, for the effective action we get 

- S eff (0; a ) « - {2) ) Y, cos (fftc - <^x + n) + O (#) , (112) 

x,n 

where 



_ A - / _ - \ / 1 — e p \ 
(2) =e-e (l-e-T) h + -^-1 • (113) 

N T 

Corrections of higher order in e P include 'abnormal' terms like cos (cp x — <£> x +ih 
similar to I 2 in ( |9~5|) . If series in a converge, corrections in a are of little 
importance in the area of N T >> 0. 



4 



1 



10 Appendix C 

We are interested in the limit N T — > oo , therefore, if there exists such 
Pasympt that for all > Pasympt then we shall have \ff > A^ o) = A min 
. Therefore, for T < T c (N T » A m i n ) we may discard all terms in ( |43] ) 
except those corresponding to j = and jo. As it can be shown a sympt — 
5/3 {gasympt ^ 1) for U(l) and I3 as ympt ^ 7/4 {g asy mpt ^ 3/2) for SU{2) 
gauge group. In all cases, which we consider, jo corresponds to fundamental 
representation, so we can preserve only two first terms in ( f43|) and may write 

- S E « A^ r In 9 (0) + exp {-iV r A G (/?)} X (0 (x)) X (« (x) f ) . (114) 
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In such a case instead of fl45| ) one may use (]35| 7.13.2(5)) 

i v ( X ) = -^=Y1 (- 2 *r m -77—^ — \ — ~ — - + o (x)- M (us) 

V2^ r (§ + v - m) r (1 + m) 

for more accurate computation of beta function, which, taking into ac- 
count ( |77D , may be written as 

Pes (f) - = -or ( ^ V 1 = - ( # In A G V 1 . (116) 



da T V dg J \dg 

In particular for G = SU (2) 



Pes (g 2 ) 




(117) 

It is easy to check that if we preserve only two first terms in ( |117| ) the beta 
function shows dip at g^ v = \/8(non-existent in this case), which, however, 
is located outside of asymptotic region. Therefore, to compute Pes (# 2 ) in- 
given approximation, we must preserve at least three first terms in (|117| ) . 
All the next terms are necessary only in the case when smooth 'link-ups' 
are needed between ( |117|) and Pes (ff 2 ) computed with the help (|116[) in a 
strong coupling area g & 1. 
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